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Abstract 

This paper presents a new, innovative procedure 
called point load superposition for determining the 
contact stresses in mating gear teeth. It is 
believed that this procedure will greatly extend 
both the range of applicability and the accuracy of 
gear contact stress analyses. 

Point load superposition is based upon funda- 
mental solutions from the theory of elasticity. It 
is an iterative numerical procedure which has dis- 
tinct advantages over the classical Hertz method, 
the finite-element method (FEM), and over existing 
applications with the boundary element method 
(BEM). Specifically, friction and sliding effects, 
which are either excluded from or difficult to 
study with the classical methods, are routinely 
handled with the new procedure. 

The paper presents the basic theory and algo- 
rithms. Several examples are presented. Results 
are consistent with those of the classical 
theories. Applications with spur gears are 
discussed . 

Introduct ion 

Gear tooth contact stress is the single most 
important factor affecting the life of a gear. 

Tooth contact stresses also have a significant 
effect upon tooth deformations and thus upon gear 
performance under load. However, in spite of their 
dominant role in gearing mechanics, contact 
stresses are still not well understood and they are 
not accurately quantified. Contact stresses are 
not nearly as well documented as, for example, root 
and fillet stresses. 

One of the reasons that gear tooth contact 
stresses are not better understood is that the 
geometry and kinematics of meshing gear teeth do 
not satisfy the assumptions required for the appli- 
cation of the classical and standard Hertzion con- 
tact stress analysis. 1 ' 2 ' 3 The Hertz analysis 
assumes a contact between frictionless cylinders. 

It does not account for noncircular geometry such 
as Involute geometry, and it does not include the 
important friction forces which occur in gear con- 
tact away from the pitch point. 

Standard finite element procedures are also 
unsatisfactory for studying contact stresses. 
Contact stresses are in essence stress concentra- 
tions which are locally nonlinear. Hence, the 


standard finite element procedures need to be 
modified through the use of gap elements, large 
numbers of elements, and iterative procedures 
before they can be accurately applied in the study 
of contact stresses. Such modifications, although 
feasible, 4 are burdensome and not generally 
available. 

A relatively new and promising procedure for 
contact stress analysis is the boundary element 
method. The boundary element method ls particu- 
larly useful In regions with high stress grad- 
ients - as with contact stresses. However, this 
method is still under development and applications 
with contact stresses also require iterative proce- 
dures to obtain the contact area geometry. 5 

In this paper we present a procedure dedicated 
strictly to contact stress analysis with direct 
application in gear tooth contact stress analysis. 
It is similar to the boundary element method. It 
does not require the restrictive assumptions of the 
Hertz analysis. That is, with this method, non- 
cylindrical geometry and frictional effects, as 
encountered in loaded meshing gears, are readily 
accommodated. The method is based upon the |>oint- 
load superposition method developed by Paul, '' 3 
and previous unpublished work of Sisera Jayasinghe 
at the University of Cincinnati. 

Historically, contact stresses have been one 
of the most important, most widely studied, and yet 
most difficult problems in the theory of elasti- 
city. Contact stress analysis can be traced back 
to the work of H. Hertz in 1881 1 ' 2 ' 3 who studied the 
contact deformation of smooth elastic bodies 
pressed together by forces normal to their contact- 
ing surfaces. He proposed that the contact region 
could be modeled by an elastic half space loaded 
over an elliptical area. 

References 9 to 22 describe recent attempts to 
use the finite element method and the boundary ele- 
ment method to determine contact stresses. 

The point-load superposition method presented 
in this paper has been outlined in a paper pre- 
sented at the Health Monitoring for Propulsion 
Systems Conference in Cincinnati ( 1990). 4,3 The 
method is also documented in a recent NASA Con- 
tractor Report of the authors. 24 

The balance of the paper is divided into nine 
parts with the following part summarizing some pre- 
liminary results useful in the sequel. The next 


three parts present the basic analysis. Numerical 
procedures are then discussed followed by some 
simple examples. The next part discusses appli- 
cation with spur gears. The final two parts con- 
tain a discussion and concluding remarks. 

Preliminary Considerations 

Consider two elastic bodies in contact. Let 
the contact region be divided into a mesh of ele- 
ments (or cells). Let the forces transmitted 
across each cell be represented by a distribution 
of concentrated forces acting on an elastic half 
space. By using long-established results from the 
theory of elasticity, together with the principle 
of superposition, the stress distribution and point 
displacements due to a cell loading are obtained. 
Then, by superposing the results from each cell the 
resulting stress and displacement distributions 
within the contacting bodies may be obtained. 

To develop the analysis consider an elastic 
half space with a concentrated normal load as in 
Fig. 1. 

Let P be a typical point within the half 
space, located by the polar coordinates (r,Q) as 
shown. Then the radial and tangential stresses at 
n 25 
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The displacement of a surface point are then 
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where c is a constant determined from a zero- 
displacement reference point. 

Next, consider an elastic half space with a 
distributed normal and tangential loading over a 
portion of the surface as depicted in Fig. 2. 

Then the resulting stresses and displacements 
at a typical point P may be obtained by super- 
position (integration) of the foregoing results. 
Following the analysis of Johnson, 25 the stresses 
at P are: 


2F n 

0 ■ - — cos u, 0 

nr 


( 1 ) 


- -J-b UzfO) (x-s)‘ 

7f 


+ g{s) (x-s) 3 ] /[ (x-s)‘ + z 2 ]1ds 


In rectangular coordinates, these stresses 
have the form: 


if 


ft* !t* 3 f 


(9) 


a = - 


2F 


* (x 2 +z 2 ) 2 


, <7 


2F 

7T 


z 3 2F xz 2 

°zx * " 

2 ... 2,2 n , .2 . _ 2 V 2 


(x +y ) 


(x +z ) 


( 2 ) 


+ z g ( s ) (x -s ) ] / { ( x -s ) 2 + z ^ ] 2 }d s 


°XZ ~ - -J! b ([Z 2 f(3)(X-S) 


+ zg(s)(x-s) j/[{x-s)‘ + ] “ )d e 


(7) 


Using the stress-strain equations and the 
strain-displacement equations, the displacements at 
a point P on the surface may be expressed as: 25 
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where E and V are the elastic constant and 
Poisson's ratio, and where h is the distance to a 
reference point of zero displacement along the 
z-axis . 


In like manner the tangential and normal dis- 
placements of a surface point are: 31 
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In like manner, if there is a concentrated 
tangential force G, the stresses at P are: 
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where c i and c are constants to be determined 
from the displacement of an arbitrary reference 
point. 


where 6 is now measured from the line of action 
of force G. In rectangular coordinates these 
stresses are 
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Application: Triangular Shape 

Loading Distribution 

Equations (7) and {8) may be used to determine 
the effects of the force distribution in the con- 
tact region. If the region is divided into a mesh 
of elements (or cells), then the effects of the 
loadings on the individual cells may be super- 
imposed. Moreover, the general loading in the con- 
tact region may be represented by the superimposed 
loadings on the cells. A triangular (that is, 
piecewise linear) loading distribution over a cell 
forms a convenient basis for modelling the general 
loading distribution. 


where C i and C 2 are constants to be determined 
from the displacements of a reference point and 
where A and B are defined as: 
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The contribution of a triangular load on an 
individual cell to the stress at a typical point 
P may be obtained by letting f(x) and g(x) of 
Eqs. (7) and (8) have the forms: 

f (x) - f 0 <a- I X I )/a , g(x) (9) 

- g 0 (a- | x | )/a , | x | $ a 


Modeling of Contact Force Distribution 


To illustrate how these results can be 
superimposed in the contact region, suppose that 
the contact force distribution is represented by 
the piecewise linear distribution, as in Fig. 4. 


where a is the half-width of the cell whose 
center is at the origin, and where f Q and g Q 
are the peak triangular normal and tangential load- 
ings as depicted in Fig. 3. 

By substituting from Eqs. (9) into (7) and by 
performing the indicated integrations, the result- 
ing stresses at P (Fig. 3) are: 


Let the contact region be divided into n 
equal-width elements (or cells) as shown. Then, 
just as the region is discretized so also are the 
force and displacement distributions. Let a be 
the element half-width. Then the element nodal 
coordinates are in multiples of a. That is, the 
x-coordinate of node k is: 
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Let u i and represent the tangential 

and normal displacements of node i due to a tri- 
angular normal load centered at node j. Then from 
Eq. (12), u and u iz may be expressed as 
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Similarly, the displacements of a surface 
point are: 
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where k is defined as: 
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Equation (15) may be used to obtain the dis- 
placements at node i due to a set of superposed 
triangular load distributions centered at each of 
the nodes. In this case the terms C ^ L and 
d if f^ are interpreted as sums over j. 
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Finally, if there is a horizontal or 
tangential loading, an analogous analysis leads to 
displacements at node i as: 
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where the coefficients c and d^ are the same 
as those of Eq. (16) and where there is a sum on j 
from 1 to n. 

The displacements of Eqs. (15) and (18) are 
measured relating to a convenient reference point, 
say the origin of the coordinate axes. 


by the friction coefficient so that alternatively 
sliding occurs. 

Figure 7 depicts the cylinders in contact with 
tangential loading. As before P^ and P„ are 
matching surface points at the initiation of con- 
tact. When the tangential load G is applied the 
distant points and Q z have rigid body dis- 
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and P 2 may be defined as: 
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Contact Analysis 

Equations (15) and (18) may be used to develop 
an iterative procedure for obtaining a detailed and 
comprehensive analysis of contact stresses and 
deformations. 

To illustrate this consider two cylinders 
pressed together as in Fig. 5. As the cylinders 
are brought together they initially have line con- 
tact. Then as they are pressed tighter the line 
contact develops into a contact strip. 

Figure 6 depicts the contact region after a 
finite deformation. The solid curves represent the 
deformed surface profiles, and the dashed curves 
represent the undeformed, or original profiles. 

Consider two points Q x and Q 2 relatively 
distant from the contact region. Then during the 
compression q and Q 2 move toward the contact 
region and toward each other through the displace- 
ments d and d respectively. 

iz 2z * 1 


If P 1 and P 2 are in the nonslip (stick) 
region is zero, and then 


(23) 


where D is the global tangential displacement of 
the cylinders. 


Numerical Procedures 

Equations (15) and (18), together with the 
constraint equations of the foregoing part, may be 
used to determine the nodal normal and tangential 
forces f and g . A difficulty which arises 
however, Is that the extent of the contact region 
is not known a-priori. Hence, an iterative pro- 
cedure needs to be developed which will determine 
the contact region as well as the nodal forces. 

For the development of this procedure, it is 
convenient to superimpose Eqs. (15) and (18) and 
express them in the form: 


Consider two typical matching surface points 
P j and P 2 separated by a distance h(x). Then 
as the deformation proceeds, P^ and P 2 simulta- 
neously move toward each other and deform inward 
into their respective cylinders. Let U and 
(measured as positive into the respective cylin- 
ders) represent the z-axis displacements of P 
and P 2 toward Q and Q 2< Then if and P 2 

come into contact, we have 

♦ Oj, ♦ h(x) - d u ♦ d,, - D, (19) 

where 0 x is the global compressive displacement 
of the cylinders. If P and P 2 are outside the 
contact region we have 

U 1M * U u * h{x) > D, { 20 > 

Equations (19) and (20) can be used to deter- 
mine the extent of the contact region. 

If there is no relative horizontal displace- 
ment of and P , we have the no-slip condition 


When the cylinders are subjected to tangential 
loading in addition to the normal loading, the 
deformation geometry depends upon the relative mag- 
nitudes of the normal and tangential loads. The 
tangential loading causes an asymmetric deformation 
pattern. This asymmetry, or distortion is limited 
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with E , V , E n , and being the elastic con- 

stants of Bodies 1 and 2. 

Similarly, for those nodes in the stick region 
Eqs. (24) and (25) together with Eq. (21) lead to 

j ~ * 0 (29) 

i*l 3*0 , . . . , ±n 


axial length and no tangential loading. The cylin- 
ders had radii: 100 and 150 in; elastic modulus: 

30xl0 6 lb/ in 2 ; and Poisson ratio: 0.3. The cylin- 
der surfaces were taken to be frictionless so that 
a comparison could be made with the classical Hertz 
solution, which assumes frictionless contact. 

Table 1 shows a comparison of numerical 
results and the Hertz solution. 
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Finally, for those nodes in the slip region we 


have 


q i =±fif i i =* -n, . . . , p, p+m, . . . , n 


(31) 


where the sign is chosen so that the force direc- 
tion is opposite to the direction of slip. 

In Eqs. (27), (29), and (31) there are 2n 1 

unknown f and 2n+l unknown g 1 (i - 0,...,±n) — 

that is 4n + 2 equation. Equations (27), (29), 

and (31) constitute 4n + 2 equations. However, 
the system is not complete since when i is zero, 
Eqs. (27) and (29) are degenerate. Hence, two 
additional equations are needed. These are 
obtained from global force summation leading to: 

F = f i A i and G = g i A i 


(32) 


where A a is the element area associated with node 
i . 


Next to examine the effects of friction (and 
hence, the consequences of neglecting friction), we 
considered two unlubricated cylinders with radii 
100 and 150 in. made of steel and aluminum respec- 
tively. Different materials were chosen to illus- 
trate the effects of material properties on the 
solutions. (The Hertz solution is based upon iden- 
tical material properties of the contacting 
bodies.) As in the first example, the cylinders 
were pressed together with a normal force of 1000 
lb per unit axial length. The elastic moduli for 
steel and aluminum were taken as: 30x10 and 

12xl0 6 lb/in. 2 and the Poisson ratios as: 0.3 and 

0.33 respectively. The friction coefficient was 
assigned as: 0.5 to simulate the unlubricated 

surfaces. Table 2 shows a comparison of numerical 
results and the Hertz solution. Figure 8 depicts 
the loading distribution in the contact region. 

Finally, to examine the effects of sliding, 
the steel/aluminum cylinders of the previous exam- 
ple were also loaded tangentially to produce 
sliding. Table 3 shows a comparison of the 
numerical results and the Hertzian solution. 

Figure 9 shows the normal pressure distri- 
bution obtained from the Hertz and numerical solu- 
tions. (Observe the shift in the numerical 
distribution due to the sliding.) Figures 10 and 
11 show the shear stress distributions on the steel 
and aluminum cylinders for the Hertz, friction and 
sliding cases. 

Application with Involute Spur Gears 


The numerical procedure is then: 

(1) Assume a contact region 

(2) Divide the contact region 
into elements or cells 

(3) Assume a stick region within 
the contact region 

(4) Solve Eqs. (27), (29), and 
(31) for the nodal forces 

(5) Adjust the extent of the 
contact region by deleting nodes 
with negative normal forces 

(6) Adjust the stick region using 
Eq. (31) 

(7) Repeat steps (1) to (6) until 
convergence is obtained 

This procedure is dependent upon superposing 
the forces from the nodes in the contact region. 
Therefore, the method is called point load 
superposition. 

Example Solutions 


To study contact stresses in mating gear teeth 
we model the teeth as contacting cylinders as 
depicted in Figure 12. The radii of the cylinders 
depends upon the position of the contact point 
along the Line of contact. 

In Fig. 12 and C n are the centers of the 

meshing gears and P is the pitch point. 1^, I 2 , 
O , 0 , P, and O are points along the pressure 
line which is tangent to the base circles of the 
gears as shown. 1^ and 1^ are at the inter- 
sections of the pressure line and lines perpen- 
dicular to the pressure line passing through 
and C respectively. and C> 2 are at the 

intersections of the pressure line and the addenda 
of the gears. O is a typical contact point along 
the pressure line. 

Then the base, pitch, and addenda circle radii 
are given by 

r bl * \^Jl\ ■ r b2 * ITO ' r pi ‘ 1^1 ' 


To illustrate the efficacy of the method we r p2 . | C 2 P | , r al - I^Oj , r a2 • | c 2 °r I 

first considered two frictionless, steel cylinders 
pressed together with a force of 1000 lb per unit 
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4 . Contact at Other Points 


where the subscripts 1 and 2 identify the associ- 
ated gear, and the notation ) C x I~ | for example, 

designates the length of the line segment between 
C and I 1 - 

Since these radii are known for a given pair 
of meshing gears, it is convenient to express the 
tooth curvatures (and hence, the contacting cylin- 
der curvatures) in terms of these radii. From the 
properties of the involute profile, the tooth radii 
of curvature at the contact point are simply the 
distances from the contact point to Z 1 and I 7 
respectively. 

Consider the following contact positions: 

1 . Contact at the Pitch Point 


In this case the radii of curvature of the 
mating tooth surfaces are: 
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where the subscripts 1 
associated gear. 

. and 2 

identify the 


2. Contact on the 

Addendum Circle of Gear 


In this case the 

radii of curvature are: 
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where <p is the pressure angle. 

3. Contact on the Addendum Circle of Gear 2 
In this case the radii of curvature are: 

P 2 - I Wl - t I Wl 2 - | C7T7I 2 ) 1/2 

(38) 

2 2 
" \ r a2 " r b2 

and 

Pi - I *7*1 1 - I W I (39) 
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Cases 2 and 3 are at the contact extremes with 
Case 1 being intermediate. For some other contact 
point, the radii of curvature will be different 
than these cases, but it can be calculated using 
the same procedures as in the above cases. 

The developed numerical procedures may now 
be used to calculate the contact stresses on the 
mating teeth. To illustrate this application, con- 
sider two identical spur gears in mesh with the 
following properties. 


Number of teeth 

18 


Pitch diameter 

3.5433 

in . 

Face width 

0.3937 

in . 

Pressure angle 

20* 


Transmitted load 

808 lb 


( tangent ial ) 
Elastic modulus 

30xl0 6 

psi 

Poisson ratio 

0.3 



The calculated maximum contact stress and 
maximum shear (Von Mises) stress for two different 
contact points are presented in Table 4. The high 
contact stress at the surface is primarily hydro- 
static loading and does not normally produce fail- 
ure. The maximum sheer stress which occurs beneath 
the surface can initiate crack formation which 
leads to pitting fatigue failure. 

Discussion 


The method employs a double discretation : one 

of the contact region, the other of the force 
distribution. The execution of the method involves 
an iterative procedure to determine the extent of 
the contact region. In this way the method is 
clearly more laborious than the classical Hertz 
method and even elementary finite-element meth- 
ods (FEM). In exchange, however, the method is 
believed to provide a more accurate and comprehen- 
sive analysis than either the Hertz or FEM 
techniques . 

Although these examples employ bodies with \ 

simple geometries (cylinders), the method outlined 
herein is not restricted to simple geometries, { 

Indeed, there are no restrictions on the geometry | 

so long as the contacting surfaces are continuous 5 

and geometrically smooth (that is, with continuous 
derivatives). j 

i 

The examples demonstrate the efficacy of the 
method. They show that the numerical results are ’ ” 

extremely close to those of the classical Hertz 
solutions, for those cases where the Hertz solution 
is applicable. Indeed, the stresses are virtually 

identical whereas there is a slight difference in ; 

the geometric results - that is, the contact region i 

width and the location of the maximum shear stress. 

(The Hertz solution assumes frictionless static 
contact between identical materials.) 
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The second and third example demonstrate the 
significant effects of friction and sliding in 
contact analysis. The examples alBo show that the 
softer material has lesser stress. 

In addition it is seen that friction and 
sliding significantly increase the shear stresses. 
Indeed, when there is sliding with a relatively 
high friction coefficient, the maximum shear stress 
occurs on the surface as opposed to being beneath 
the surface (Figs. 10 and 11). 

The conditions of gear tooth contact (Non— 
cylindrical surface, friction, and sliding) do not 
satisfy the assumptions for Hertz theory. The 
point load superposition method is not restricted 
by these assumptions. The gear stress results of 
Table 4 are believed to be accurate, limited only 
by the precision of the descretat ion . 

The method outlined herein is developed for a 
two-dimensional analysis with a triangular load on 
the individual elements or cells. The triangular 
load provides for a piecewise linear representation 
of the contact loading. The method can be extended 
to higher order load representations (for example, 
the use of cubic splines) and to three-dimensional 
analysis . 
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A point load superposition procedure based on 
the theory of elasticity has been used to study 
contact stresses between elastic bodies, including 
meshing spur gear teeth. This method is validated 
by comparison with results from the classical Hertz 
method for frictionless cylinders in contact. The 
following specific conclusions are reached: 

1. A new method of contact stress analysis 
has been presented. It is a numerical method which 
has fewer restrictions than classical procedures. 

2. The method is applicable for a broad range 
of geometries of the contacting bodies. The con- 
tacting bodies may have different material pro- 
perties. The contacting surfaces need not be 
smooth. Sliding between the contacting surfaces 
is also permitted. 

3 . The accuracy of the method is demonstrated 
by a comparison with results from classical methods 
for simple cases where the classical method is 
applicable. 

4. The established efficacy of the method 
justifies applications with contacting gear teeth. 
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Table 1. - Comparison of Hertz and Numerical Solutions for 


Contacting Steel Cylinders without Friction 


Item 

Hertz 

Solution 

Numerical 

Solution 

Difference 

Width of Contact Region 
(in, ) 

0.136156 

0.137287 

0.83% 

Maximum Normal Pressure 
(psi) 


9358.90 

0.08% 

Maximum Shear Stress (psi) 

2808.05 

2810.26 

0.08% 

Location of Maximum Shear 
on Z-Axis (in.) 

0.05310 

0.05352 

0.79% 


Table 2. - Comparison of Hertz Solutions (Identical, Frictionless 
Materials) and Numerical Solutions for Contacting 
Steel/Aluminum Cylinders with Friction (The numbers 
in parenthesis refer to the aluminum cylinder) 


Item 

Hertz 

Solution 

Numerical 

Solution 

Difference 

Width of Contact Region 
( in. ) 

0.18896 

0.18873 

0.12% 

Maximum Normal Pressure 
(psi) 

6738.12 

6810.89 

1.08% 

Maximum Shear Stress (psi) 

2021.43 

2081.32 

(2031.28) 

2.95% 
(0.48% ) 

Location of Maximum Shear 
on Z-Axis ( in . ) 

0.07369 

0.07967 

(0.06461) 

8.1% 

(13%) 


Table 3. - Comparison of Hertz Solutions (No Sliding) and Numerical Solu- 
tions for Contacting/Sliding, Steel/Aluminum Cylinders 


Item 

Hertz 

Solution 

Numerical 

Solution 

Difference 

Width of Contact Region (in.) 

0.18896 

0.19041 

0.77% 

Maximum Normal Pressure (psi) 

6738.12 

6737.87 

0.00% 

Maximum Shear Stress (psi) 

2021.43 

3365.39 

66.4% 

Location of Maximum Shear 
on Z-Axis ( in. ) 

0.07369 

0.00 

100% 

Shift of Center (in.) 

0.00 

0.01123 

100% 


Table 4. - Gear Tooth Stress Results and Comparisons 



Pitch Point 
Contact 

Addendum Contact 

Surface Contact Stress (psi) 

194 , 480 

286,013 

Maximum Shear Stress (psi) 

58 , 407 

85,899 

Depth to Maximium Shear 
Stress (in.) 

0.0056 

0.0038 
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Figure 7. Contact Region with Tangential Loading 



a => half of contact width Trom numerical solution 
P 0 = maximum normal pressure from numerical solution 

normal pressure distribution 

a tangential pressure distribution 

(directions are indicated by arrows) 


Figure 8. Contact Region Loading Distribution for 
Contacting SteeL/Aluminum Cylinders With Friction 
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a - half of contact width from Hertzian contact 
R, = maximum norma! pressure from Hertzian contact 
a sliding contact 

Hertzian contuct 


Figure 9. Comparison of Hertz and Numerical Solutions 
for Normal Pressure for Contacting/Sliding, 
Steel/Aluminum Cylinders 



r maximum shear stress 
P Q maximum normal pressure (Hertz) 
a half of contact width (Hertz) 


Figure 10. Shear Stress Distributions 
on the Steel Cylinder. 
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T maximum shear stress 
P 0 maximum normal pressure (Hertz) 
a half of contact width (Hertz) 


Figure 11, Shear Stress Distributions 
on the Aluminum Cylinder 



Figure 12. Modeling Mating Gear Teeth as Contacting Cylinders 
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